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Nomenclature 

0 Anglo of attack, deg 

S cs Control surface deflections, deg (+down) 
f Emissivity 
p Density, slug/ ft 5 

<t Stefan-Boltzmann constant, 4.76 x 10 ft 2 R l — .s 
b Wingspan, 27.7 ft 
h Altitude, kft 
L Overall length, 58.3 ft 
M Mach number 
p Pressure, lb/ft 2 
q Heating rate, Btu/ft 2 -s 
Rt Reynolds Number 
T Temperature, °F 

1 Velocity, ft /s 

,r Fuselage station, ft 
y Buttline station, ft 
z Waterline station, ft 


Introduction 

As early as May 1993, Orbital Sciences Corporation (OSC) of Dulles, Virginia began 
considering low cost methods to launch relatively small payloads (~1,500 lbs) into space. 1 ' 2 In 
the fall of 1993, OSC approached NASA to inquire about an alliance which would utilize the 
Agency’s space transportation knowledge base, facilities, and analysis tools. Subsequently. 
Presidential Decision Directive NSTC-4 led NASA to announce a Cooperative Agreement 
Notice 1 to build a small reusable launch vehicle (R.LV) designated X-34. 1 ’ 2 ' 1 fi This was 
awarded to OSC, teamed with Rockwell International, in the Spring of 1995. By February 
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199G, however. OSC and Rockwell had determined that venture would not be profitable 
and withdrew from the agreement. 7,8 In the Spring of 1996, NASA solicited proposals for a 
different, vehicle, also designated X-34. OSC was awarded this contract in June of 1996. 

The new X-34 is to be a suborbital technology demonstrator, capable of flying to Mach 
8, reaching an altitude of 250,000 ft, and landing autonomously on a conventional runway. A 
typical mission profile is shown in Fig. 1. The vehicle is launched from OSC's L- 1011, ignites 
its liquid rocket engine, ascends to altitude, and then (-oasts to a down-range landing/' The 
vehicle is to serve as a test bed for a multitude of RLV technologies such as composite airframe 
and propellant tank components, low-cost avionics via the Global Positioning System (GPS), 
and a flush air data system. 1011 

This paper presents results used to define and design the Thermal Protection System 
(TPS) necessary for the vehicle's entry into the Earth's atmosphere. This work is part of a 
larger effort used to design the TPS, including engineering codes, arc-jet facilities, and wind- 
tunnel studies. 12 1 ' The TPS is predominately insulation blankets 1, of the Flexible Reusable 
Surface 1 Insulation (FRSI) class with SIRCA tiles 10 used only on the leading edges. 

NASA Langley's aerothermal contributions consisted of several focused efforts. This 
paper describes three 1 ’ benchmark-quality, full- vehicle, flight condition computations using 
state-of-the-art. Computational Fluid Dynamics (CFD) to ascertain TPS temperature limits. 
Riley et a/. 12 anchored to these two results, expanded the trajectory envelope using a coupled 
inviscid-boundary layer method. Incorporating the data from this limited array of trajectory 
points, Wurster et al 11 developed a time history of heating rates for the entire trajectory 
using engineering codes. This analysis yields the heat load for over sixty points on the 

rt A very similar mission profile to the X-15 nearly forty years previous. 9 

h An additional case which does not include the aft portion of the vehicle is also presented. 
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vehicle, providing data necessary to determine the appropriate TPS thickness. In addition 
to the above. Berry at al lA created an extensive database with several wind tunnel entries, 
investigating effects of changing Mach number, configuration, angle of attack. Reynolds 
number, and control surface deflections. 

The remainder of this paper begins with a discussion of the numerical tool and geometry 
employed, and is followed by a brief discussion of the trajectory and the selected points. 
Next are the results, consisting of: a detailed description of the dominant flow features, a 
comparison with experiment in the form of surface' shear patterns, and surface temperat ures, 
heating rates, and pressures. Finally, a discussion of TPS design issues are discussed, followed 
by some concluding remarks. 


Numerical Method 

In the hypersonic continuum regime, the Langley Aerothermodynamic Upwind Relax- 
ation Algorithm (Laura) CFD tool 18,19 is used to describe the aerothermodynamics of X-34. 
Laura is an upwind-biased, point-implicit relaxation algorithm for obtaining the numeri- 
cal solution to the Reynolds-averaged Navier-Stokes equations for three-dimensional viscous 
hypersonic flows in thermo-chemical nonequilibrium. 20 The upwind-biased inviscid flux is 
constructed using Roe’s flux-difference-splitting 21 and Harten’s entropy fix 22 with second- 
order corrections based on Yee’s symmetric total-variation-diminishing scheme. 22 This is the 
same computational tool that has been used to describe? the aerothermodynamics of blunt 
body shapes such as Mars Pathfinder, 24 25 Mars Microprobe, 26 Stardust, 27 and COMET; 28 
and other vehicles such as X-33, 29 HL-20, 20 Space Shuttle, 21 22 and Reentry-F. 24 

For all of the results contained within, the Laura code was run assuming air to behave as 
a perfect gas and the full Navier-Stokes equations were slightly simplified via the thin-layer 
assumption (see Ref. 35 for rationale). Furthermore, for nearly all the results presented, 
the flow was assumed to be fully turbulent, using the Balwin-Lomax algebraic turbulence 
model 2 '’’ 1 ' modified with a damping term according to Gupta et a/ 2(> as implemented by 
Cheatwood. 2, However, selected laminar results an* also presented to bound the problem. 

The wall-temperature boundary condition is specified as the radiation-equilibrium, wall 
temperature according to the Stefan-Boltzmann relation, 

q = mT\ (1) 

The radiative-equilibrium wall temperature, T, is coupled during the solution procedure with 

( A higher-order turbulence model (e.#., Spalart-Allmaras or k-jj) which may be more appropriate for the 
massively separated flow on the lee side of the vehicle, was not available in Laura at the time of this study. 


3 



tlio wall heating rate, </, where f is the surface emissivity' 1 and a is the Stefen-Bolt/mann 
constant. 

This coupled solution procedure between wall temperature and wall heating rate yields 
accurate temperature and heat ing rates because of the good insulation features of the selected 
TPS blankets. 


Geometry 

The X-34 geometrical description used in this study is designated X0001215 as received 
from OSC in Initial Graphics Exchange Specification (IGES) format. Figure 2 shows the 
X-34 geometry as modeled. Overall, the vehicle is similar to the Space Shuttle, having a 
cranked delta wing planform and vertical tail. Note, however, that the fuselage transitions 
to rectangular cross section approximately midway back and that the wing terminates before 
reaching the aft end of the vehicle. The vehicle is to weigh 45,000 lbs with full fuel and have 
a wingspan, b. of 27.7 ft and an overall length, L, of 58.3 ft. 

The coordinate axes are defined in the typical body-oriented manner: x running longitu- 
dinally. a along the starboard wing, and z defined by the right-hand rule, pointing upward 
with the origin located at the nose. 1 ' 

Small modifications were made to the original geometry to allow more tractable grid 
generation and facilitate obtaining the CFD solutions. These modifications consisted of the 
following: 

• Backward or forward facing steps created by TPS material interfaces such as that 
caused by the wing leading edge SIRCA tiles and the trailing AFRSI blankets were 
not preserved as sharp steps. The surface grid lines were mapped onto the stepped 
surface, but no attempt was made to align them with the discontinuities. As a result, 
the sharp steps art* replaced by ramps. 

• The span-wise gaps between the elevons and the body and the elevens and the ailerons 
were filled in as solid surfaces. Any other gaps or surface irregularities have not been 
modeled. Past experience has shown that when elevons and ailerons are slaved to move 
together as is the case for the hypersonic portion of the X-3 Fs flight profile, that the 
effects of the gap between them is a highly localized phenomenon as shown by Berry 
r.t. al . u However, filling the gap between the body and the elevon has a larger effect 
because it channels flow onto the aft portion of the vehicle which would normally have 
passed through to the lee-side wake. 

‘'A constant emissivity of 0.8 was used for all surfaces. 

‘OSC uses a different origin location. 
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b) Wind side view. 

Figure 2. X-34 geometry. 
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Figure 3. Stagnation-point heating (hot-wall) for three X-34 trajectories. 

• The lower, aft portion of the tail surface overhanging the hack of the vehicle was clipped 
off to provide' a continuous supersonic outflow boundary condition. 

• The vehicle' wake and base re'gion including the engine' nozzle' and lee-side of the body 
flap we're neither modeled nor computed. 

The' surface and volume meshes were generated using Icemcfd,** 8 GridTool, 49 Grid- 
GEN. ,() 3DMAGGS, 41 and Vgm 12 as described by Alter. 14 A typical full-vehicle grid including 
the body flap and wing wake has a total of 70 blocks and 9 million grid points. 

Trajectory Cases 

Figure' 3 shows stagnation-point heating rates, hot- wall, for a 7.11 inch radius sph(Te f as 
predicted by Fay and Riddell theory for three' different X-34 trajectories provided by OSC. 
The' XI 004701 trajectory was used by Palmer et a/ 15 and Milos et a l u> to pre'dict the heating 
environment to be' encountered by the nose cap and wing leading edge SIRCA tiles. The 
X1004601 trajectory was used as the reference heating trajectory for the overall TPS design 
per OSC's request. Tin' actual flight trajectory is anticipated to be closer to XI 003904. 

Table 1 presents the free-stream conditions and control surface deflections for the eases 
computed. R Tin' points shown in the Table are taken from the X1004601 trajectory and 

correspond to peak nose cap heating and slightly thereafter when the angle of attack is 
significantly lowered. 

f The nose radius of the X0001215 configuration. 

R Note: for Cast 1 2h the solution was not computed aft of the wing trailing edge. 
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Table 1. Free-stream conditions and co ntrol surface deflections 


Case 

tics 

(deg) 

h 

(kft) 

Moo 

(-) 

a 

(deg) 

Voc 

(ft/s) 

Po 0 

(slug/ft 3 ) 

Too 

(»F) 

Re oo 

(-) 

la 

lb 

0 

10 

118.4 

6.3 

23 

6490 

1.42xl0' 5 

-22 

16x10® 

2a 

2b 

0 

10 

112.1 

6 

15.2 

6110 

1.89xl0^ 5 

-30 

21xl0 6 


Results 

First, the dominant flow field features are presented and analyzed to provide a basis for 
the discussion of the results to follow. Next, a comparison with experimental data of Berry 
et a/ 13 is given for Case 2 conditions in the form of near-surface streamlines via the oil-flow 
technique. Surface distributions of temperature are presented next, followed by a sample of 
the associated heating rates and an example of the effects of a laminar or turbulent boundary 

layer. 

Note: discontinuities in the figures are a result of block-to-block averaging errors when 
a block-marching solution strategy was used. These errors are aesthetic only, the compu- 
tational cell-centered values are globally conservative and consistent across block-to-block 
boundaries. 

Dominant Flow-field Features 

Figure 4 shows four stream-wise cuts with contours of pressure non-dimensionalized by 
free-stream, p, for Case lb. Figure 4(a) also includes red lines indicating the location of the 
other three stream-wise cuts and, in blue, the locations of span- wise cuts that will be shown 
in Figs. 5 and 6. 

As shown in Fig. 4(a), the forebody of the vehicle is enveloped in a highly swept bow 
shock which intersects with the wing-induced shock. As the flow passes the wing tip the 
ou tor wing shock swoops aft ward duo to the expansion. 

Figure 4(b) depicts a stream-wise cut along the fuselage outboard of the centerline. The 
wind-side bow shock is shown clearly. Also shown are the faint remnants of the shock created 
by the deflected elevons; and, at the aft end of the vehicle, the embedded shock due to the 
deflected body flap. The lee-side shock is also apparent, although somewhat smeared due to 
a highly stretched grid in this region. Toward the front, on the lee-side, a small expansion- 
recompression is evident due to a TPS transition. This is followed by a rapid expansion over 
the canopy region. A second, slight compression is also noticeable due to cross flow effects 
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b) Vertically, along fuselage (2y/b=0.15). 


Figure 4. Non-dimensional pressure contours for various stream-wise cutting planes 

(Case lb). 
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d) Vertically, along wing (2y/b=0.60). 


Figure 4. Concluded. 
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from the side of the fuselage. This is followed by further expansion, until the point where the 
fuselage transitions from a rounded cross-sectional shape to one with a considerably smaller 
corner radius (see Fig. 2(a)). At this juncture, a shock is generated, which is then followed 
by a small expansion when the transition is complete. The wake of the wing is responsible 
for the area of lowest pressure just ahead of the tail which again raises the pressure slightly 
due to an embedded shock. 

Figure 4(c) shows a stream-wise cut further away from the centerline, but still well inside 
of the wing-bow shock interaction. The shock generated by the deflected control surface is 
readily apparent on the wind side. Another dominant feature present is the wake of the wing 
and strake which is altered by a fish-tail shock structure emanating from the trailing edge 
of the wing. Moving outside of the bow-shock interaction area as shown in Fig. 4(d) shows 
similar features. 

Span- wise cuts for Case lb are shown in Figs. 5 and 6. h Figure 5 shows cross-flow 
streamlines colored by pressure' and Fig. 6 shows Mach number contours from 0.2 to 9.0 in 
increments of 0.2. 

The first pair of sub-figures (5(a) and 6(a)) are taken on the forebody, just behind the nose 
cap. A pair of wind-side vortices is clearly seen near the centerline. There is no evidence of 
vortices on the lee side at this axial station. The over-expansion and subsequent compression 
is evident along the side of the fuselage. 

A very different picture is evident in Fig. 5(i) which shows a cross-section located near 
the wing trailing edge. There is no longer any evidence of the wind-side vortices, and a pair 
of well-defined vortices sit on the lee side of the fuselage. Small wing-tip vortices are present, 
and the wing-strake vortices lie along the side of the fuselage. The source structure above 
the wing at mid-span is merely the demarcation of inboard/outboard span-wise flow as well 
as the location where the flow has either a positive or negative vertical component. Also 
apparent is a small separation bubble on the inboard wing lee side. 

On the wind side the embedded shock due to the deflected control surface is also clearly 
indicated by a merging of streamlines paralleling the deflected control surfaces. Figure 6(i) 
also shows the presence of the extensive hour-glass shaped cross-flow shock standing on the 
fuselage. This shock serves to turn the flow parallel to the vehicle symmetry plane. The 
vertical gradients in Mach number between the cross-flow shocks (which are not present in 
pressure see Fig. 5(i)) are believed to mark the entropy interface between flow that has gone 
through only the bow shock and that which has gone through both the bow and cross-flow 
shocks. 

h Recall, the relative vehicle locations are depicted by the blue lines in Fig. 4(a). 

'No explicit scale is given; as with Fig. 4, blue indicates low pressures while red indicates high pressures. 
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i) Before wing trailing edge (x/L=0.79). 


j) Behind wing trailing edge (x/L— 0.86). 


Figure 5. Concluded. 


Next, the flow field in the vicinity of the body flap is examined with Figs. 7 and 8. 
Looking from the underside of vehicle, Fig. l shows the view in the vicinity of the body flap 
notch region. The vehicle extends off to the right and the rest of the body flap continues 
to the left. The back surface of the vehicle is truncated due to the limited computational 
domain used to model this region. Near-surface streamlines are shown colored by pressure 
variation. A re-attachment line is clearly' evident along the leading edge of the body flap 
notch. This is created by the flow expanding around the bottom, aft corner of the vehicle. 
In fact, as shown by the series of span-wise cross-sectional cuts of density contours in Fig. 8, 
the flow is strong enough to create an embedded shock along the leading edge of the body 
flap notch. For this figure, the viewpoint is from the top of the wing, looking aft towards the 
side of the body flap (the body and body flap surfaces have been removed). The first cross- 
section contains the aft end of the vehicle, with the following cross-sections in the notched 
area, and the last cross section just after the body flap regains its full span. 

Note: since the grid topology in the notch region was such that it was not possible to 
run the notch block with the algebraic turbulence model as implemented, this block was 
run using the laminar Navier-Stokes equations. As a consequence, any surface wetted by 
flow emanating from the notch region has erroneous surface temperatures. However, from 
investigations comparing the purely laminar results and the mixed turbulent-laminar results, 
the basic flow structure in this area remains un-effected. This is fortunate since Fig. 7 shows 
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a) Downstream of spherical nose cap b) Canopy region (x/L=0.11). 

(x/L=0.04). 




c) Before strake begins (x/L=0.29). 


d) On strake (x/L=0.39). 


Figure 6. Mach contours (AM=0.2, Case lb). 
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Figure 7. Streamlines colored by pressure in the vicinity of the body flap notch 

(Case la). 
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Figure 8. Span-wise cross-sectional cuts showing density contours along the aft of the 

vehicle and the body flap (Case la). 

that there is a dividing streamline between the flow which remains on the wind side and the 
flow which travels onto the side wall and then returns to the wind side. This serves to limit 
the contamination caused by the embedded laminar block. 

Surface Shear Patterns — Comparison with Experiment 

Computer-simulated near-surface streamlines at flight conditions are shown in Fig. 9 
for Case 2a, assuming a turbulent boundary layer. Figure 10 shows the corresponding oil 
flow from a Mach 10 wind-tunnel run where the attached flow over model is assumed to 
be laminar. 13 While Mach 6 wind tunnel results are also available, the Mach 10 data was 
chosen for comparison since the normal-shock density ratio is more comparable to the flight 
condition. 

Comparing wind-side flow features, both datasets show similar indications of wind-side 
vortex patterns (refer to Fig. 5(a)) scrubbing the boundary layer toward the centerline on 
the forward quarter of the vehicle. Forced by a span-wise pressure gradient generated by 
the nose geometry which transitions from a spherical cross section to a rounded, triangular 
shape, the wind-side flow is initially converging toward the centerline. After a quarter of the 
vehicle, the centerline boundary layer has thickened and these wind-side vortices dissipate, 
establishing the streamline patterns containing an outboard component that dominate the 
surface flow over the latter three-fourths of the vehicle. The surface discontinuity at the 
wing-fuselage junction gives rise to a slight change in the streamline patterns, more visible 


17 



a) Wind side. 



b) Starboard side. 



c) Lee side. 

Figure 9. Computed near-surface streamlines (Case 2a). 
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in th(' experimental results. Both datasets show a divergence of streamlines emanating from 
the wing-st rake joint which is next to the how-shock impingement region. Tin' experimental 
streamlines at the wing trailing edge show some localized behavior due to the elevon gaps, 
not present in the computational results, which did not model the elevon gaps. 

Side-view comparisons show very good agreement, particularly in the forebody patterns. 
Both datasets show a strong cross-flow shock location and the separation line on the fuselage 
< a used In the vortex emanating from the strake. A strong interaction region is seen, partic- 
ularly in tin' computational solution, on the fuselage aft. of the wing trailing edge. There is a 
bleed-through effect in this region on the experimental model due to the elevon gap. which 
is not modeled in the numerical simulation. 

Lee-side surface features are in good agreement between the datasets, showing the strong 
cross-flow shock separation on the fuselage and the strake vortex across the wing, close to 
the fuselage. Inboard flow is seen over much of the wing while the fuselage shows outboard 
How. generated by a longitudinal recirculation zone. Turning and separation, with a probable 
horse-shoe vortex, occurs in front of the vertical tail. 

Further comparisons with experimental data are available in Ref. 13. 

Surface Temperatures 

This section presents the surface temperatures for all turbulent boundary-layer cases 
computed: Cases la. lb, 2a, and 2b. Discussion includes the effects of angle-of-attaek 
variation and control surface deflections. 

Case la: Figure 11 shows surface temperature contours' for Case la. The wind-side 

view shown in Fig. 11(a) shows the expected high temperatures for the nose stagnation 
region. The forebode chines also show significant temperature elevation due to the vehicle’s 
forebody geometry approaching the bow shock generated by the spherical nose cap. Also 
present on Fig. 11(a) are streaks of elevated temperatures crossing the wings chord-wise. 
This corresponds to the remnants of the forebodv bow shock (cf. Fig. 4(a)). Also, although 
not readily discernible with this choice of temperature contour levels, there is a second streak 
emanating from the bow shock-wing shock interaction point just outside of the wing/st rake 
junctme, inked aft from tin? leading edge of the wing. This second streak is readily apparent 
in the experimental results of Ref. 13. The vast acreage of the wind-side fuselage surface is 
in the range of 1300-1400 °F. 

Note that the outside edges of the body flap have slightly cooler temperatures than the 
inboard portion. This area of cooler temperatures is associated with near-surface streamlines 

J Note that there are three contour-level legends in the figure, one for each of the lee-side, starboard, and 
wind-side views; but these same three contour level distributions are held constant for all subsequent surface 
temperature figures to facilitate comparisons. 
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that emanate from the leading edge of the notched body flap. This is trend is erroneous as 
discussed on page 13. 

The starboard and lee-side views (Figs. 11(b) and 11(c)) show the effects of several flow 
features previously discussed. For instance, between the canopy and the round-to-squared 
fuselage transition, longitudinal low-temperature streaks are well defined, corresponding to 
the pair of vortices located just above the fuselage. As mentioned previously on page 10, 
when the fuselage cross-section transition is encountered, an embedded shock is generated 
" is indicated by the triangular-shaped region of higher temperatures in this area. 

Case lb: Surface temperatures for the deflected control surface case (lb), are shown in 

Fig. 12. As compared to the undeflected case, the temperatures on the forebody are not 
effected bv the change due to fact that the shock layer has supersonic flow. As shown bv 

Fig. 12(a) the deflected surfaces now experience 200-300 °F higher temperatures on their 
wind sides. 

For this case, the non-physical effect of having to run the “notched” body flap region 
with the laminar equations is more readily apparent than for the undeflected case since the 
embedded body-flap-notch leading-edge shock (refer to Fig. 8) is stronger. 

Although difficult to see in Figs. 11(b) and 12(b), the wing wake is slightly larger for the 
deflected case. Also, the footprint of the wing fish-tail shock is at a lower angle relative to 
the vehicle’s longitudinal axis. 

Case 2a: Surface temperatures for this lower-angle-of-attack case are shown in Fig. 13. 

As compared with Case 1, the wind-side temperatures are on the order of 250 °F lower 
while maintaining the same qualitative distribution, except for the wing-shock/bow-shock 
interaction occurring slightly closer to the fuselage. Figs. 13(c) and 13(b), however, show 
increased heating on the lee-side. Most notable are the nose region and the tail leading edge. 
The fuselage transition area is also slightly hotter. The wing-wake structure appears to have 
higher energy as its temperature footprints are on the order of a 100 °F higher than the 
higher-angle-of-attack case. The terminating fish-tail shock also appears stronger, reflected 
by its more oblique angle. 

Case 2b: The effects of deflecting the elevon control surface 10 deg. are shown in Fig. 14. 

As for Cases la and lb, the forward portion of the vehicle is not effected by this change due 
to supersonic flow, and the deflection simply results in higher wind-side temperatures for the 
flap (see Fig. 14(a)). However, for this lower angle-of-attack case, there is a small area of 
flow separation near the outboard end of tin' elevon, resulting in slightly lower temperatures 
in this area. This is faintly visible in Fig. 14(a) as a sliver of blue along the hinge line 
between 70 and 80% span. Its presence is more clearly indicated by the convex shape of the 
yellow-green contour level on the elevens in this region. 
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T (°F) 


200 400 600 800 1000 1200 1400 



b) Starboard side. 




c) Lee side. 


Figure 12. Surface temperatures (Case lb). 
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a) Wind side. 



b) Starboard side. 



c) Lee side. 

Figure 13. Surface temperatures (Case 2a). 
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T (°F) 



a) Wind side. 


T (°F) 



b) Starboard side. 



c) Lee side. 


Figure 14. Surface temperatures (Case 2b). 
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Heating Rates 

Figs. 11-14 have been presented in terms of temperatures since this is the most important 
quantity to consider when ensuring that no TPS temperature limits are violated. For the 
reader, estimates of the corresponding heating rates can be derived from the temperature 
distributions with Eq. 1. 

An example of an actual heating rate distribution is shown in Fig. 15 k for Case 2b. The 
heating-rate distributions for the remaining cases are presented in Figs. 16-18. 

Boundary-Layer State 

For Cases la and lb, both laminar and turbulent boundary-layer solutions were obtained. 
This section provides a brief synopsis of the effects of the boundary-layer state, focusing on 
the deflected control surfaces case. Case lb. 

Figures 15 and 19 show surface heating rates for turbulent and laminar boundary layers, 
respectively. By comparing Figs. 15 and 19, it is immediately evident that the turbulent 
boundary-layer assumption provides a more conservative estimate of the heating rates than 
the laminar boundary layer. Overall, the laminar heating rates tend to be half of those for 
the turbulent boundary layer. However, the qualitative agreement is very close, signifying 

that most of the heating patterns are determined by “inviscid” flow features such as shocks 
and vortex cores. 

Laminar heating results for Case la are presented in Fig. 20 while Figs. 21 and 22 show 
laminar temperature distributions for Cases la and lb. 

Figures 23 26 show turbulent pressure distributions for Cases la, lb, 2a, and 2b, respec- 
tively while Figs. 27 and 28 show the laminar results for Cases la and lb. Upon comparing 
Figs. 24(a) and 28(a) for Case lb, only small changes in the pressure distributions on the 
wind side are evident. These are primarily due to the fact that the laminar boundary layer 
sustains a small separated region in front of the deflected control surfaces whereas the turbu- 
lent boundary layer does not. This moves the associated embedded shocks forward. On the 
lee side. Figs. 24(b), 24(c), 28(b), and 28(c) show the major differences to be in the footprint 
of the wing-strake vortex structure and the lack of the wing fish-tail shock signature on the 
side of the fuselage for the laminar case. 

Tiie dashed pattern of very low heating rates apparent on the lee side of the fuselage (see Fig. 15(c)) 
ts due to an instability of the Baldwin-Lomax turbulence model when coupled with the iterative, radiative- 
equilibrium wall boundary condition in a multiple-block context. An improved algorithm for determining the 
location of the maximum vorticity within the boundary layer eliminated this problem in subsequent runs. 
This particular case was not re-computed due to limited resources and the fact that most regions of interest 
are outside the domain of dependence. 
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a) Wind side. 



b) Starboard side. 



c) Lee side. 

Figure 15. Heating rates (Case lb). 
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a) Wind side. 



c) Lee side. 

Figure 16. Heating rates (Case la). 
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a) Wind side. 



b) Starboard side. 



c) Lee side. 


Figure 17. Heating rates (Case 2a). 
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a) Wind side. 


q (Btu/tf-g) 

0.1 0.3 0.5 0.8 1.1 1.4 3.2 7.9 10.0 



b) Starboard side. 



c) Lee side. 


Figure 18. Heating rates (Case 2b). 
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a) Wind side. 



b) Starboard side. 



c) Lee side. 

Figure 19. Heating rates (Case lb), laminar. 
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a) Wind side. 



b) Starboard side. 



c) Lee side. 


Figure 20. Heating rates (Case la), laminar . 
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a) Wind side. 



b) Starboard side. 



c) Lee side. 

Figure 21. Temperature Distributions (Case la), laminar . 
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a) Wind side. 



b) Starboard side. 



c) Lee side. 

Figure 22. Temperature Distributions (Case lb), IcLTninar. 
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a) Wind side. 



b) Starboard side. 



c) Lee side. 

Figure 23. Pressure Distributions (Case la). 
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a) Wind side. 



b) Starboard side. 



c) Lee side. 


Figure 24. Pressure Distributions (Case lb). 
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a) Wind side. 



b) Starboard side. 



c) Lee side. 

Figure 25. Pressure Distributions (Case 2a). 
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a) Wind side. 
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b) Starboard side. 



c) Lee side. 

Figure 26. Pressure Distributions (Case 2b). 
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a) Wind side. 



b) Starboard side. 


> 



c) Lee side. 


Figure 27. Pressure Distributions (Case la), laminar. 
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a) Wind side. 



b) Starboard side. 



c) Lee side. 


Figure 28. Pressure Distributions (Case lb), laminar . 
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Figure 29. Thermal protection system multi-use temperature limits (°F). 


TPS Considerations 

As provided by OSG, the general layout of the TPS materials is given in Fig. 29. indicating 
the materials’ respective temperature limits while the TPS blanket orientations are depicted 
in Fig. 30. 

As a result of comparing the Case 2 temperature distributions from Figs. 13 and 14 and 
the TPS temperature limits shown in Fig. 29, the bond line between the AFRSI and FRSI 
blankets on the lee side near the wing landing edge was moved aft to account for the increased 
exposure found at 15 dogs, angle of attack. The TPS treatment on the lee side of the nose 
region was also changed for this reason. 

By comparing predicted near- surface streamlines of Fig. 9, the TPS material layout in 
Fig. 29, and the TPS blanket orientations of Fig. 30, one can investigate aspects of the 
TPS orientations. For example, upon carefully examining these figures, it is evident that no 
blanket-to-blanket or blanket-to-tile gaps are aligned with the local flow. 


Concluding Remarks 

Radiative equilibrium surface temperatures, heating rates, streamlines, surface pressures, 
and flow-field features as predicted by the Laura code were presented for the X-34 Tech- 
nology Demonstrator. Results for two trajectory points near peak heating and the effects of 
control surface deflections were presented for fully turbulent flow. Laminar flow results were 
also presented to illustrate the effects of the boundary-layer state. The effects of the results 
on the TPS design were also indicated. The presented results of wall temperature were used 
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Figure 30. Thermal protection blanket layout. 

extensively to ensure that the temperature limits of the selected TPS blankets will not be 
exceeded. 

Results show that downward deflection of control surfaces substantially increases the 
temperatures they experience during flight. Furthermore, a “notched” body flap design 
creates an embedded shock which creates large heating rates on the tapered portion of the 
body-flap side wall. 

The effect of lowering the angle of attack from 23 to 15 deg. at flight conditions was shown 
to substantially increase the extent of lee-side heating on the nose-canopv region and aft of 
the wing leading edges. Also shown was that using turbulent boundary-layer heating results 
is considerably more conservative than using results which model a fully laminar boundary 
layer. 
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